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1 Introduction 

Modern application of quantum mechanics increased interest in the space of quantum states: positive operators normalized 
by the assumption that their trace is fixed, Tr p — 1 . For applications in the theory of quantum information processing it 
is often sufficient to restrict the attention to the operators acting on a finite dimensional Hilbert space. 

The set Mn of density matrices of size N forms a convex body embedded in K/'" -1 . In other words it forms a cross- 
section of the cone of positive operators with a hyperplane corresponding to the normalization condition. In the simplest 
case of one qubit the set Mi is equivalent, with respect to the Hilbert-Schmidt (Euclidean) geometry, to a three dimen- 
sional ball, B3. For higher dimensions the geometry of Mn gets more complicated and differs from the ball B AT 2_ 1 [1,2]. 

Non trivial properties of the set of mixed quantum states attracted recently a lot of attention. The volume V, the 
hyperarea A and the radius R of the maximal ball inscribed into Mn was computed with respect to the Hilbert-Schmidt 
measure [3], which leads to the Euclidean geometry. The volume of the set of quantum states was computed with respect 
to the Bures measure related to quantum distinguishability [4], and also with respect to a wide class of measures induced 
by monotone Riemanian metrics [5, 6]. The set Mn is known to be of a constant height [7], so the ratio A/V coincides 
with the dimensionality of this set, equal to N 2 - 1 . 

If N is a composite number, the density operators from Mn can represent states of a composed physical system. In this 
case one defines the set of separable states, which forms a convex subset of the set of all states, A1^, ep c M N . Although 
a lot of work has been done to estimate the volume of the subset of separable states [8, 9, 10, 1 1, 12, 13], the problem of 
finding the exact value of the ratio Vol (vV(^ ep ) / Vol (vVl^) remains open even in the simplest case of two qubits [14]. 

In parallel with investigation of the set of quantum states, one studies properties of the set of completely positive (CP) 
maps which act on Mn- Such maps are important not only from the theoretical point of view: for instance linear CP 
maps acting on a two-level quantum system correspond to linear optical devices used in polarisation optics [15]. Due to 
Jamiolkowski isomorphism [16, 17], the set of trace preserving CP maps acting on Mn forms a N 4 - N 2 dimensional 
subset of the N 4 - 1 dimensional set Mn 1 - of states acting on an extended Hilbert space, < Hn® c ^n- In the simplest case of 
N — 2 the structure of this 12-dimensional convex set of maps was studied in [18]. 

We start this paper by reviewing the properties of the set of subnormalized states for which Tr p < 1 . Such states are 
obtained by taking the convex hull of the set of normalized states and the particular "zero state". In the classical case, one 
could argue that such a step is equivalent to increasing the number of distinguishable events by one, and renaming into 
N + 1 . This reasoning is based on the fact that the set of subnormalized states of size N as well as the set of normalized 
states of size N + 1 form TV-dimensional simplices. However, this is not the case in the quantum set-up, in which the 
set of subnormalized states M s ^ b has TV 2 dimensions, in contrast to TV 2 + 27V dimensional set Mn+i- The fact that the 
dimensionality of the set of subnormalized states grows with the number V of distinguishable states exactly as V 2 plays a 
key role in an axiomatic approach to quantum mechanics of Hardy [19]. 
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The main aim of this work is to describe the set of completely positive trace non-increasing maps which act on the set 
Mn of mixed states. We compute the exact volume of this A^ 4 dimensional set with respect to the Euclidean (Hilbert- 
Schmidt) measure and characterize its extremal points. The trace non-increasing maps have a realistic physical motivation, 
since they describe an experiment, for which with a certain probability the apparatus does not work. This could be an 
interpretation of the "zero map" after action of which no result is recorded. Such maps are sometimes called trace 
decreasing [20], but to emphasize that the set of these maps contains also all trace preserving maps, we prefer to use a 
more precise name of trace non-increasing maps, (TNI). 

The paper is organized as follows. In Section 2 we analyze the set of subnormalized states and compute its volume. The 
measures in the set of mixed states induced by partial trace are investigated in Section 3. In Section 4 we define the set 
of trace non-increasing maps and provide its characterization, while in Section 5 the volume of this set is calculated with 
respect to the flat (Hilbert-Schmidt) measure. The set of extremal trace non-increasing maps is studied in Section 6. 

2 Subnormalized quantum states 

Let Mn denote the set of normalized quantum states acting on jV-dimensional Hilbert space 

Mn '■= {p ■ p 1 = p, p> 0, Trp = 1} • 

Mn forms a convex set of dimensionality (n 2 - l). In the simplest case N = 2, this set is equivalent to the Block ball, 
Mi = B 3 c R 3 . 



M^ := [a- : o 4 = cr, cr > 0, Tr cr < 1} = conv hull {0, M N ) ■ (1) 



Definition 1 

An Hermitian, positive operator cr is called a subnormalized state, if Tr cr < 1 . 

The set of subnormalized states acting on the jV-dimensional Hilbert space 'Hn will be denoted by M?^° . By construction 
this set has N 2 dimensions and can be defined as a convex hull of the zero operator and the set of quantum states (see 
Fig. 1), 

For instance, the set M?^ forms a four-dimensional cone with apex at and base formed by the Bloch ball Mi — B3 . Note 
that the dimension of A1™ b grows exactly as squared dimension of the Hilbert space. Due to this fact the subnormalized 
states are a convenient notion to be used in an axiomatic approach to quantum theory [19]. 

In order to characterize the set Al^ ub of subnormalized states we compute in this Section its volume with respect to the 
flat H S-measure induced by the Hilbert-Schmidt metric [3] . Consider the set Mi of 2 x 2 density matrices, parameterized 
by the real Bloch coherence vector £ € B3, 

I 2 -> -> 

P=y+^-H , (2) 

where S denotes the vector of three rescaled traceless Pauli matrices Note that with such a 

normalization the radius of the Bloch ball B3 is given by Ri = 1 / V2, as it can be obtained from the relation Tr p 2 < 1 . 
With this definition, the HS-distance between any two density operators, defined as the HS (Frobenius) norm of their 
difference, proves to be equal to the Euclidean distance De(£i,£2) between the labeling Bloch vectors £1,^2 £ B3 c R , 



D Hsfe.Pg) = ^Tr[( P| _ - P| J 



fl-fe = D E(fl,ft) • (3) 



The above formula holds for an arbitrary N, provided that 



p=^+^-£ , (4) 
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Figure 1: The set of subnormalized states: (a) the set of eigenvalues for N - 3, (b) the convex cone of A^ 2 dimensions 
with zero state at the apex and the (n 2 - l) dimensional set Mn as the base. 



the real coherence vector £ in (4) is taken (n 2 - l)-dimensional and H now represents an operator-valued vector which 

consists of (n 2 - l) traceless Hermitian generators of SU (AO, fulfilling Tr (H,E 7 ) = d (J . Note however that in this case the 

geometry of the space of coherence vectors £ does not coincide with the ball B A/ 2_ 1 , but constitutes instead a convex subset 
of it [2]. The condition Trp 2 < 1 yields the upper bound for the length of the coherence vector, |£| < ^(N _ 1)/N -■ Rn- 
In the case N - 3 the vector E consists of the set of 8 normalized Gell-Mann matrices {S,}f =1 . 

A metric space consisting of a set Mn and a distance d is automatically endowed with a measure induced by the metric: 
The measure is defined by the assumption that all balls of a fixed radius defined in Mn with respect to the distance d have 
the same volume. 

The infinitesimal HS-measure around any matrix p e Mn factorizes as [3] 



d /"HS (P) = TTT dv (Ai, . . . , Ajv) X dv 



Haar 



AH 



(5) 



where the factor dv A represents the measure in the simplex A^_i of eigenvalues Ai, . . . , A N , while dv Haar depends on the 
eigenvectors of p. The pre-factor V^V emerges in (5) when we force the variables Ai, . . . , A/v to live on the simplex \n-i ■ 
This reduces the number of independent variables by one, sinceJl^! A, = 1, and introduces a factor -y/detg in (5), where 
g denotes the metric tensor in the (Af - l)-dimensional simplex. Such a metric arises due to a change from Af linearly 
dependent variables {A,}^ to the {N - I) linearly independent ones {A,}^ 1 . 

The last factor dv Haar can be integrated on the entire complex flag manifold [21, 3] defined by the coset space, Fl^ :- 
U(A r )/[U(l)]' v , that is the space of equivalence classes of matrices U diagonalising the given p. The volume of the flag 
manifold induced by the parametrization used in (4) is given by [3] 



VolK] = X 



dv 



Haar 



0/2 



«8° 



1 ! 2! ■ ■ ■ (AT — 1) ! 



(6) 
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Even after splitting off the A-phases of [11(1)]", a residual arbitrariness still remains in the diagonalization of p, related 
to the fact that different permutations of N generically different eigenvalues A, belong to the same unitary orbit. This 
explains the factor N ! in (5), given by the number of equivalent Weyl chambers of the simplex Ajv-i . 
The measure dv A (Ai , . . . , A N ) reads [3] 



I N 

dv A (A u ...,A N ) = d 



( N \ N _ 

2 a- - 1 n ® (Ai) n ( A; - a >) 2 dAi • ■ • aan ■ (?) 

w=l / 1=1 i<j 

where the Dirac delta and the product of the Heaviside step functions ensure that the measure is concentrated on the 
simplex An-i- Up to a normalization constant, the r.h.s. of the above equation defines a probability distribution on the 
simplex, 

dAAi,...,A w )oci»®(Ai Atf)dAi ...dA N , (8) 

where the upper index 2 refers to the exponent in the last factor of equation (7). In random matrix theory it is commonly 
denoted by ft and called repulsion exponent, as it determines the repulsion between adjacent eigenvalues. It is equal to 1, 
2 and 4 for real, respectively complex, respectively symplectic ensembles of random matrices [22]. In this work we are 
going to restrict our attention to complex density matrices so we fix /3 = 2, but one may also repeat our analysis for other 
universality classes. 

For later purpose, we now introduce a bigger family of probability distributions, indexed by a real parameter a: 



P ( f\A u ...,A N ):= C { f ] 6 



1-£a, f]0(A,)Af- 1 f](A,-A ; ) 2 , (9a) 



i=l / i'=l i<j 

with normalization constant C^' 2) given by 

^ N n£i m + j) nj +(<*-!)] 



i r ( N \ N 2 

—— : = I 5 1-VA/ r~f ©(A,) Af _1 I""! (a, - A ,) dA x . . . dA N - 



(9b) 



The probability distribution in (8) represents a special case of P ( N ' , being 

P ( V(A U ...,A N ) = P ( f\A u . . ., A N )\ a=l . (10) 
Putting together equations (5-7) and (9b) one derives [3] the (HS-) volume of the set Mn of mixed quantum states 



As shown later in this Section, the volume of the set Al^ ub of subnormalized states can be easily obtained using the (flat) 
Euclidean geometry as the volume of the cone with base given by Mn ■ However, for later use, we shall first derive this 
formula directly by integrating an extension of the HS-measure over the entire cone. Let us first formulate the following 
Lemma, proved in Appendix A 

Lemma 1 

Consider a one parameter family of probability measures dv#(Ai, . . . , A N ) defined on the set CH (N) := 
conv hull {0, \n-\] 

N _ 

dv K (Ai ,...,A N ) = Y] ©(A,) Af- N Y] (A, - Ajf dAi . . . dA,v (12) 
i=l i<j 
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and labeled by integers K > N. Then the volume v K (CH (AT)) reads 

v*(CH(A)) = f dv J ,(A 1 ,...,A JV ) 
Jo 



;ch(ao KNCf- N+l » ' 

where C^~ N+l ' 2 ^ is the coefficient defined in (9b), with a = K - N + 1. 
The measure d/^^Js' ^ on tne set °^ surj normalized states Al^ ub is given by 

dy« S Hs(^) = ^7dv sub (A 1 ,...,A w )xdv Haar , (13) 

and differs from the one of equation (5), relative to quantum normalized states, by the factor yfR. In equation (5), such a 
factor was due to the change of variables needed to express the volume elements in terms of the set of (N - 1) independent 
variables on the simplex Ajv-i; in the present case we do not need to change the variables anymore, and the factor VA 
does not appear. 

Moreover, the definition of the measure dv sub on the entire set CH (A) of Lemma 1 used in equation (13) differs from 
the one defined on the simplex Ajy_i and used in equation (7). In particular, fovK-N equation (12) implies 



dv sub (Ai, ...,A N ) = dv N (A u ...,A N ) = ® 



N \ N _ 

^a,-i n ®( a ') n ( A - - a j) 2 dAi • ■ • aan ■ ( u) 

i=i ) i=\ i<j 



In analogy to the derivation of the volume (1 1) of the set of normalized states, the HS-volume of the set Al^, ub of subnor- 
malized states can be earned by equations (13, 14, 6) and (9b) together with Lemma 1, in which we set K = N: 

Vol hs (AC) = >< Vol K1 >< >* < CH W> = W NM2 ^^TiN^ ' (15) 

The above result can be obtained directly using the Archimedean formula for the Euclidean volume of the D-dimensional 
cone of Figure 1 (panel b) 

V=^A-h , (16) 

where V is the D-dimensional volume of the (hyper-) cone representing Al^ ub , A is the area of its (D - l)-dimensional 
base Mn, and h denotes its height, that is the distance between the base and the apex (the latter corresponding to the state 
(To = 0). Making use of the definition of HS-distance (3), one gets the results 

V =Vol H s(A^ ub ) 

D = Dimension (A^, ub ) = N 2 

A =Vo1hs(A1 w ) • ( 17 > 

h = D HS (My, 0) = inf D HS (p, 0) = 

peM N -yTV 

Note that h = Dhs (p*, 0), where p+ = Ijf/N denotes the maximally mixed state of Fig. I, due to the following chain of 
relations: 

N 1 
h 2 = inf Z)L (p, 0) = inf Trp 2 = inf Y A 2 = inf |a 2 + ■ ■ ■ A 2 J = — ■ (18) 

peM N HS ^' peM N H AeA„_, ^ A]+ -+A W =1 L 1 N \ N 

Substituting results of (17) into (16) we arrive at 

( i,sub\ VoIhs (Mat) 
Vol H s(Al w )= Wl (19) 

which, due to (1 1), is consistent with the volume given by (15). 
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2.1 Generating random subnormalized mixed states with respect to HS measure 

As we have already seen, the HS-measure endows Mn with the flat, Euclidean geometry. Moreover, as emphasized in 
FIG. 1, the set A1^ ub of subnormalized states constitutes with respect to this geometry an jV 2 -dimensional cone, whose 
(N 2 - l)-dimensional base is Mn- Thus, directly from (1), every state <x e At^, ub can be decomposed as 

<r = ap + (l-a)0 , (20) 

where p e Mn, is the null-state and a is a positive number less or equal to 1 . Therefore generating states {<x, } uniformly 
in the cone means to distribute homogeneously p, in Mn and to use them in (20). The weights a, e [0, 1] must be 
distributed accordingly to a density function f(a) scaling as a N _1 . These random numbers may be obtained by inverting 
the cumulative distribution function F(x) := J f(x)dx over uniformly distributed random numbers e [0, 1], that is 

generating homogeneously e [0, 1] and taking a,- := F _1 (£) = %V N . As a result, given a sequence of HS-distributed 
mixed states {p,} e Mn , as the one studied in Section 3.1, and a sequence {£,-} of random numbers uniformly distributed 
in [0, 1], one obtains an algorithmic prescription to generate a sequence of HS-distributed random subnormalized states 

3 Measure induced by partial trace of mixed states 

The aim of this Section is to determine the measure induced on a bipartite K x N quantum system AB, represented by 
means of an extended Hilbert space Hab = Ha ® Hb = C K ® <C N , by partial tracing over one of the two subsystem. Such 
measures will be essential in determining the volumes of various sets of maps, that is done in the subsequent Sections 
of this work. Without loss of generality, we consider ancillary systems A whose dimension K is greater or equal to the 
dimension N of the system B. The induced measure depends on the choice of the systems to trace over, and on the initial 
distribution of the density operators p AB on the composite system "Hab ■ We start by analyzing a propaedeutic example: 

3.1 Partial tracing over pure states 

Consider random pure states pab = \^ab){4>ab\ distributed according to the natural Fubini-Study (FS) measure. Such 
a measure is the only unitarily invariant one on the space of pure states, and for N = 2 it gives the measure uniformly 
covering the entire boundary of the Bloch ball [2]. A general pure state can be represented in a product basis as 

k n 

\Vab) = X X Mji ^ A ® ^ B ' 

''=1 7=1 

The positive matrix p B = MM^ is equal to the density matrix obtained by a partial trace on the Zf-dimensional space A, 
Pb = Tr A (pab)- The spectrum of p B coincides with the set of Schmidt coefficients A,- of the pure state \i//ab)- The matrix 
M needs not to be Hermitian, the only constraint is the trace condition, Tr MM* = 1, that makes p B a density matrix. 
Furthermore, the natural measure on the space of pure states corresponds to taking M from the Ginibre ensemble [22] and 
then renormalize them to ensure Tr MM 1 = 1 [23]. The probability distribution of the Schmidt coefficients implied by the 
FS measure on Ha ® 1~Ib is given by [24] 

i - Yj a <- n ® (Ai) k ^ n n ( Ai _ Aj ) 2 ' (2ia) 

i ) i i< j 



P% K (A U ...,A N ) = B% K 6 
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in which the upper index 2 is the repulsion exponent, as we deal with complex density matrices. The normalization 
constant B^ K reads [23] 

d(2) TjKN) 

N ' K '~WJ^r(K-j)T(N- j + i) ' {nb) 

Observe that the entire distribution (9), including the normalization constant, can be tuned into eq. (21), provided we 
choose a - 1 = K - N. In particular, for a = 1, that is K = N, equation (10) shows that the induced distribution for 
the eigenvalues of pB coincides with the HS-distribution. Thus, generating normalized Wishart matrices MM 1 ', with M 
belonging to the Ginibre ensemble of N x N matrices, becomes a useful procedure for producing N xN density matrices 
HS-distributed, and the algorithm described in Section 2.1 becomes effective. 



3.2 Partial tracing over random HS-distributed mixed states 

Consider now a related problem of determining the probability distribution of states p B = Tr A (pab), as it is assumed that 
the mixed states pab are distributed according to the HS-measure on Mkn- 

From Section 3.1 we know that p A B itself can be generated by taking FS-distributed pure state Pa'B'ab = Wa'B'Ab){^A'B'Ab\ 
in an extended Hilbert space "Ha'B'ab = *Ha'B' ® "T^ab = & KN ® C KN and by partial tracing over the ancillary subsystem 
A'B'. Putting all together we obtain 

p B = Tr 4 \pab] = Tr A [Tr A ^ (Pa'B'ab)] = Ttaa'B' [\&a>b>ab)(&a>b'Ab\] ■ (22) 

The latter equation implies that the desired distribution is obtained by coupling the A-dimensional system B to an en- 
vironment of size NK 2 , generating FS-distributed pure states in this overall N 2 K 2 system, and finally tracing out the 
environment. Hence the distribution of the spectrum of ps is given by P ^(Ai, . . . , A#) of equation (21). 

In the rest of this Section we will consider the special case K = N. This assumption corresponds to generating mixed 
states distributed according to the HS measure on the space At// 2 of bipartite systems, and tracing out one of them. In 
particular we will focus on the simplest cases of N = 2 (qubits) and N = 3 (qutrits). 



3.3 Partial trace of two-qubits mixed quantum states 

Let us start with the simplest case N = 2, for which B^\^ = B®\ = 180 180 and the simplex Ai, corresponding to the 
positive Ai, A2 such that A] + A2 = 1, is nothing but an interval [0, 1]. The probability distribution (21) reads in this case 

f^(Ai, A 2 ) = 180 180 6(1 - A, - A 2 ) 0(A t ) ®(A 2 ) A* A^ (A, - A 2 ) 2 • (23) 
We express the eigenvalues Ai and A 2 in terms of a real parameter r e \-\, 5] = : &i as follows 



1 

A 2 = I-r 



and we earn from (23) the radial distribution inside the Bloch ball 



1 \ 6 

1 .21 „2 



where %■? (r) denotes the indicator function of the simplex ts. x (see FIG. 2a). 



(24) 



P(r) = 720720 {r)\--T*\r* , (25) 
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3.4 Partial trace of two-qutrits mixed quantum states 

For N = 3 the eigenvalues Ai, A2 and A3 can be expressed in polar coordinates (r, (f) in a form similar to (24), 

A] = i + rcos(^ + l^-J 
• A 2 = \ + rcos(0) • (26) 

A3 = I + r cos {cf> - ^tt) 

Similarly as before, we indicate with A2 the counter-image of the simplex A2, that is the set in the (r, <p) plane such that 
(Ai(r,0,A2(r,0,A 3 (r,0))eA2. 

Computing the Jacobian of the transformation (26) (with A3 = 1 - Ai - A2) we see that the volume element transforms as 

V3 

dV = dAi dA 2 = — x X2 2 (r,</>) x rdrd</> . 

The value of the radial variable r is related to the purity of the mixed state, 3/2 r 2 — Tr p 2 - 1 /3, where Tr p 2 = A 2 + A 2 + A 2 . 
The constant B^ N3 follows from (21b) 

B(2) = 80! 
3,27 12 -24! -25! -26! 
We are now in the position to compute for this case the explicit probability distribution (21) 

Pf 21 (A u A 2 , A 3 ) = B% ^(Ai, A 2 , A 3 ) A 24 A 24 A 24 (A 1 - A 2 ) 2 (Ai - A3) 2 (A 2 - A3) 2 , (27) 

that in the polar plane reads 

— 9 80' /ll 1 \ 24 

P(r ' 0) = 64 24! -25l -26! ^ ^ * °® (27 " 4 ^ + 4 ? ™ " ^ 

The latter distribution is invariant under the transformations (/>—>-(/> and </>—></> + 2kn/3, k e Z, as shown in FIG. 2. 



4 Trace-non-increasing maps 

A generic state of a TV-dimensional quantum system is completely described once a positive, Hermitian and normalized 
density matrix p e Mn is given. In order to analyze and classify the set of all possible physical operations on a quantum 
system, we need to describe the set of superoperators <b mapping Mn onto itself. The linearity of physical operations can 
be expressed by forcing the superoperator <D : Mn Mn to act as a matrix action on the "vector" p, that is 

p' = Op or p' mfl = <& mfl p„ v 

nv 

We use Einstein summation convention for indices appearing twice. In order to map the domain Mn onto itself, the linear 
super-operator <J) must fulfill these additional 

Properties 1 

i) P' = (P') f 

ii) Trp' = Trp =1 

iii) p > 



<=> O m// = <D*„ m (29a) 

nv vn 

«* <l>Mm = 6 nv (29b) 

nv 

<^ Qmiipnv > when p > • (29c) 
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Figure 2: Induced distributions in the simplices of eigenvalues obtained by partial trace of mixed states ofNxN systems 
distributed according to HS measure, (a) N = 2, A\ = [-1/2,1/2]. Comparison of an histogram with theoretical 
prediction (25) represented by solid line, (b) N = 3. Contour lines of the distribution (28) in the simplex A2. (c) N = 3. 
Plot obtained numerically for 10 5 random states of size N 2 -9. 



In particular, all the superoperators <£> fulfilling property (29b) are called trace preserving (TP). For any given <5, it proves 
convenient to introduce the dynamical matrix D®, uniquely and linearly obtained from the superoperator by reshuffling 
the indices [25, 26]: 

(30) 



D mn = <D 



' mn 





<=> 


D mn — D 

fiv mn 


SO 




(31a) 


Trp' =Trp =1 


<=> 


D mn — & nv 

mv 


SO 




(31b) 


p' >o 


<=> 


D mn P nv 
iiv 


when 


P>0 


(31c) 



The dynamical matrix can be thought as a matrix on a bipartite N xN quantum system AB, represented by means of an 
extended Hilbert space 'Hab of the same kind of the one in Section 3, so that D mn = A {m\ ® B {n\ \fi) A ® |v) fi . 

tiv 

In terms of Do, Properties 1 can be re-expressed as follows: 
Properties 2 

i) 
ii) 
iii) 

In the following we will focus on the set CPn of maps O : ^(C^) h-> ^(C^) that are completely positive [27], that is on 
maps O such that for every identity map I K : C K h-> <D K the extended maps O ® I K : S(C N ) kC^h ^(C^) ® are 
positive (here £(<C N ) is the Banach space of bounded linear operators on C N ). A very important property characterizes 
the dynamical matrix of a completely positive (CP) map <b: 

Theorem (Choi [28]) : A linear map <t> is completely positive if and only if the corresponding dynamical 
matrix is positive. 

Note that property (31c) holds true in general once that the positivity of the matrix D® is proved for product states of 
<C N ® <C N , as stated in the Jamiolkowski Theorem [16]. This property is implied by complete positivity, that is a stronger 
condition. If we combine equations (31a-31b) with the complete positivity, we observe that for any given O e CPjf 
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(i.e. the set of completely positive trace preserving maps), its rescaled dynamical matrix p® := D®/N possesses the 
three properties of Hermiticity, positivity and normalization. Thus p® belongs to the set M N2 of density matrices of size 
N 2 . This leads to the celebrated Jamiolkowski isomorphism [16], that maps the entire set CP N of quantum maps onto a 
proper (N 4 - jV 2 )-dimensional subset of the set M Nl of quantum states on extended system. This subset, denoted by M 1 ^ 
, contains all states p such that Tr A p = In- The N 2 constraints reducing the dimension of M N2 come from equation (31b). 
With the aim of removing these N 2 constraints, we introduce a family of linear maps: 

Definition 2 

A linear positive map <£ is called trace-non-increasing (TNI), if Tr 0(p) < Tr p = 1 for any p e Mn ■ 

We state now a Lemma that makes a link between CP-TNI maps and their images given by Jamiolkowski isomorphism, 
that is the set of their (rescaled) dynamical matrices: 

Lemma 2 



For any given <t> e CP N , its rescaled dynamical matrix cr® := D®/N is Hermitian (according to (31a)), 
positive (as in the statement of Choi Theorem) and fulfills the following constraint: 

ii') TrO(p) = Trp' < Trp = 1 o Tr^ cr® < ^ • (32) 

For a proof of this Lemma, see Appendix B. Note that the above constraint is a kind of a relaxation of (31b). The reduced 
and rescaled dynamical matrix Tr^cr^, belongs then to a subset of the set defined in (1). To specify this subset 

consider the following set of sub-tracial states 

M D N := (cr e M% ub : cr < ^} = [cr e M s N ub : max [EV (cr)] < 1} , (33) 

where EV (cr) denotes the set of eigenvalues of cr. Note that the set M N coincides, up to rescaling by N, with the set of 
positive Hermitian operators £, fulfilling £, < . A collection of such operators, which satisfies the relation 2, = Lv , 
is called POVM (Positive Operator Valued Measure) and plays a crucial role in the theory of quantum measurement [2]. 
Definition (33) allows us to rewrite condition (32) as 

cD eCP T N ™ TvAO-^eM^ ■ (34) 

Finally, we can summarize the result of equation (34), in the next 

Proposition 

Every trace non increasing map O e CP™ can be represented by a sub-tracial state o- e M N c Al^, ub , 
whose explicit expression is given in terms of the rescaled dynamical matrix cr of Lemma 2 by 

o- := Tr A cr® = ^Tr A D® ■ 

Moreover, the rescaled dynamical matrix cr® offers itself another representation of cD into an A^ 4 -dimensional 
proper subset of Al^ u 2 b , as one can derive from the right hand side of (32) that 

Tr cr® = Tr AB cr® = Tr B (Tr A cr®) < Tr B (I N /N) = 1 ■ 

In the next Definition 3, we will denote this set with M® 2 . 
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Map O 


Rescaled dynamical 
matrix cr® 


Reduced and rescaled 
dynamical matrix <xo 


CP Trace Preserving maps 








CP Trace Non Increasing maps 


CP™ 







Table 1: The two kinds of CP-TP and CP-TNI maps analyzed in this Section are mapped in the set of their correspondent 
superoperators <E>, the set of their rescaled dynamical matrices ctq, and the set of their representative sub-tracial states &®. 
Inclusion relations between the sets under consideration are explicitly shown in the table. 



Definition 3 

M m N2 := [cr e M s N f : Tr A a e M n N } . 

To clarify the notation used we collect the sets of states and maps considered in table 1 . 

In the space of the N-eigenvalues of NxN positive Hermitian matrices, the set of sub-tracial states defined in (33) 
consists of a cube inscribed into the set Al^ ub of subnormalized states of size N. Such a multi-dimensional cube has a 
vertex in the origin, is oriented along axes, and touches the simplex of quantum states Mn in a single point I/v/W, as shown 
in FIG. 3 for N = 2 and N = 3. Observe that every cr® of M N can be rescaled in order to let it belong to the cube of At° . 
In terms of maps it means that every <D e CP N can be mixed with the map in order to become trace-non-increasing. 

The set M^ 2 of CP-TNI maps has N 4 dimensions. It contains the N 4 -N 2 dimensional set At^ of CP-TP maps, which 
satisfy Tr <D(p) = Tr (p), as a subset. The set CP T N P forms extremal points in CP™. 



5 Volume of the set of CP TNI maps 



In Section 4, the Jamiolkowski isomorphisms allowed us to establish links between superoperators <5 e CP™ 1 , rescaled 
dynamical matrices cr = Tr^ cr e At^ , and sub-tracial states <xo e Ai^. In this Section we will make use of this 
representation aiming to define a measure in CP™ and compute its volume. 

We start defining the Hilbert Schmidt measure d//^ s on the space of trace non-increasing maps O e CP N as the HS 
measure of their representative reduced dynamical matrices cr<j, e Al^ 2 c Al^f. Such a measure is analogous to the one 
introduced in (13-14), restricted from Al^ u 2 b to M® 2 , and reads 



M m cA1 sub 

N 2 N 2 



1 



dv sub (A! ,...,A^) 



N 2 N 2 



x dv 



Haar 



(35) 



Here A : , . . . , K N i denote the eigenvalues of the matrix cr® of size A^ 2 . 
The restriction of the space from Al^ u 2 b to Al^, 2 does not affect the Haar measure on the entire complex flag manifold 

Fl ( ^ \ given by the coset space U(A r2 )/[U(l)] A ' 2 of equivalence classes of unitary matrices of eigenstates of cr . The 
volume of this flag manifold follows from (6) 



Vol 



H 



<> 



dv 



Haar 



(2n) 



N 2 (N 2 -l)/2 



1 ! 2 ! ■■■(N 2 -l)l 



(36) 
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Figure 3: The set of CP trace-non-increasing maps O e OP™ represented, by means of the Jamiolkowski isomorphism, 
as the subset M° of the set of subnormalized states M™ h , in terms of the reduced and rescaled dynamical matrices 
(t,d = Tr A cr . The eigenvalues of these matrices <x<p span an jV-dimensional cube inscribed into the set of subnormalized 
states of size N, here plotted for (a) N - 2 and (b) N -3. 



In order to parallel (11) and (15) we need to compute the volume spanned by the eigenvalues of matrices contained in 
Al^ according to the measure dv sub (A! , . . . , A N i), multiply it by the volume (36) of the flag manifold Fl ( £ ' and divide 
the result by N 2 \ as in (35). Performing this task one obtains 



Vol hs {OPT) = Vol [EV (M^)] x Vol [f/f >] = JL f dv sub (A, , . . . , A N2 ) x Vol 



Fl (N) 



(37) 



Note that the operation of partial trace maps the set M® 2 into M° . Assuming that mixed states from M® 2 c Al^ u 2 aie 
distributed according to the HS measure, we need to analyze the measure induced in M° c A1^ ub by partial trace. 
In Section 3.2 we described how the normalized HS-distribution on Mnk, whose expression is equivalent to the induced 



distribution P^ KNK (A] , . . . , A^) of equation (21), is mapped by partial tracing into the distribution P^' NK 2 (Ai > • • • > A;v)- 



,(2) 



In particular, for K — N, the HS-distribution on Aiffi is mapped into P NN , (Ai, 

In a completely equivalent way, the measure dv sub (Ai, . . . , A N i) = dv^Ai, 
by partial tracing into the measure dv^Ai, . . . , Ajv), and (37) yields 



■ ,A N ). 

,A N i) of equation (35) is transformed 



Vol HS (C^ NI )= t4t f , ,dv JV 2(A 1 ,...,A w2 )x\bl 

N 1 JEV (^2) 



fi {N) 



Vol 



Ft 



A 2 ! 



I 



EV(M°) 



dv N ,(A l ,...,A N ) 



(38) 



Making use of (12) and (36) we arrive at 

VV 2 (A< 2 -l)/2 ^ N 



vo1 hs (cpt) = iTz,...^, f n a7 _a ' * n & ^ a ?~ n n ( a <- - ^ ^ ■ • • - 
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where we have made explicit the domain of integration using the Heaviside step function 0. Last integral can be computed 
using next Lemma 3, whose proof is in Appendix C. 

Lemma 3 

On the A/-dimensional cube a N := |0, ^ , consider the one parameter family of measures dv°(A l5 . . . , A;v) 



N I 1 \ 2 

dv° (Aj , . . . , Ajv) = Y\ ® [jj ~ A, 0(A,) Af- N Y\ (A, - A,) . . . dA N 

i'=l ' »"</ 



(40) 



labeled by any integer K > N. Then the volume v° (a N ) reads 



C J □ / A A N + 1,1,1, AQ .... 

J^dv K (A u ...,A N ) = — , (41) 



where I(K - N + 1, 1, 1, N) is the Selberg's integral [22] 



Finally, using Lemma 2, equation (39) yields the final formula for the HS volume of the set of TNI maps 



( 7n \N 2 {N 2 -\)l2 r 



_ (InfWjW J^ Tjl+j) r(N 3 -N + j) r(j) 

~ 1!2! ■ ■ - N 2 I H N^T^+j) (43) 

Equation (38) implies that the HS volume of the set of CP-TNI maps is proportional to the volume of the set of sub- 
tracial states At^ of (33) computed according to the induced measure dv^Ai , . . . , An). In particular Lemma 3 allows 
us to compute such a volume, which is in turn proportional to the volume spanned by all possible elements of a POVM, 
according to any given induced measure dv^(A 1 , . . . , A N ). Moreover, as we explained in the Remark of page 17, also the 
volume of all sub-normalized states A1^, ub can be computed according to the any induced measure. For comparison, we 
present below their ratio for any value of parameter K, reminding the reader that in the former calculations we set K — N 3 . 

-^M = C<™ X I(K - N + 1, 1, 1, AT) = ^ fl -I®- • 

6 Extremal CP TNI maps 

We start this Section by introducing another representation for CP-TN I maps O : p i — > p', alternative to the one described 
in Section 4, and summarized in the next Lemma 4 (for a proof see Appendix D). 

Lemma 4 

For any given CP-TNI map O, its action on density matrices p e Mn can be described by using a discrete 
family of N x N operators A^, whose number S does not exceed N 2 . The explicit action of the map O reads 

s 

<5>(p) = Y, A nP A l < 44a > 
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and the matrices fulfill 

s 

^AlA^N&l , (44b) 

where the sub-tracial state o"o e At^ is given by the reduced and rescaled dynamical matrix correspondent 
to the map as in the Proposition of pag. 10, and cr T means transpose of cr. Conversely, if S ^ N 2 operators 
Af, of size NxN fulfill 

1 S 

then 

s 

®(p):=Yj A »P A l ( 45b > 

defines a CP-TNI map. 

When O is a CP-TP map, that is when cr^ = crj = In/N, then equations (44) define the usual Kraus-Stinespring 
representation and the operators A p are said Kraus operators. In this sense Lemma 4 can be seen as an extension of the 
KS-representation. 

We aim to use the previous Lemma on a particular sub-class of CP-TNI maps, represented by very trivial sub-tracial 
states, as in the following 

Definition 4 

We say that a map O e CP™ is k-extremal if its representative sub-tracial state <xo e is diagonal and 
possesses exactly k non vanishing entries, all equal to 1/N. Those sub-tracial states represent maximally 
mixed states in A:-dimensional subspace. In the following such maps will be denoted by (b^. 

Let £l N denote the ordered set of positive integer less or equal to N, and let P (£!#) be the set "sorted parts of Qjy", that 
is the set of all possible ordered collections of elements of £l N . Consider the partition 

N 

V{£l N )= [jp k (Q N ) , (46) 

k=Q 

where each Vk (Cl N ) contains exactly k elements, Po (Q^) = {0} and P N (Q^) = \Qn)- Thus it proves convenient to use 
the natural isomorphism 

fe*» t (n*) ~ *e = jfT i \ e)(el (47) 

for labeling the set of k-extremal CP-TNI maps: see Figure 4 for an example with N = 3, on which 

P 3 (a 3 ) = {£l 3 } = {{1,2,3}} , P 2 (Q 3 ) = {{1,2U1,3},{2,3}} , P x (£1 3 ) = {{1},{2},{3}} and P (Q 3 ) = {<d} . 

The vectors \ () are meant to be orthonormal . Figure 4 illustrates another geometrical meaning of the integer number k 
labeling a A:-extremal map: it is proportional to the taxi distance in between and the vector of diagonal entries of 
Of , that is the minimal number of sides connecting to &^ in the plot of eigenvalues. 

From (47) we earn that <Xf is essentially a projector operator rescaled by N, &( = PJN: we denote as Atf the image of 
Mn obtained by projection, 

M c := |J P (P P ( <zMT . (48) 

peM N 

Using the isomorphism (47), Lemma 4 can be re-phrased as follows: 
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As 




Figure 4: Set of CP-TNI maps with extremal points marked. In particular (*) marks the CP-TP map $3 = I3/3, (❖) 
marks the map Oq = 0, (<?) and (*) marks maps <E>2, respectively Oi . 



Lemma 5 

For any given A:-extremal CP-TNI map <t>£, that is for any given £ e (Q-n), its action on density matrices 
p € Mn can be described by using a discrete family of NxN operators A^, whose number S does not exceed 
N 2 . The explicit action of the map reads 

s 

<S>(p) = Yj A »P A l (49a) 

n=\ 

and the matrices fulfill 

s 

YjAlA^N&s . (49b) 

Simply using the previous Lemma 5 we deduce the following 
Properties 3 

Let <!>£ be a ^-extremal CP-TNI map, represented by some &^ , with £ e Tk (&n)- Then it holds true: 

1) when k — N, <1>n belongs to OPjf '. When k < N the trace preserving condition becomes p-dependent: 
nevertheless <t>£ e CP N acts as a Trace Preserving map on p e M% ; 

2) acting on the maximally mixed state p* := In/N reveals the parameter k, being Tr (p*)] = k/N; 

3) A1^ ub 2 (S>k (p*) = (k/N) to, for some to e Mn- When k = N and O is unital, then to = p* . 
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Proof of Properties 3: 

1) Tracing both sides of equation (49a), using the cyclicity property of the trace, and inserting (49b), one obtains 
Tr [Ojt (p)] = Tr <x^p) = Tr (P(p), where P( is the projector of equation (48). Then using the property of projector 
and ciclicity, together with definition (48), one gets 

Tr [O k (p)] = Tr [p] p) = Tr (P ( p P ( ) . (50) 

Assuming that p e M( means, according to definition (48), that it exists a quantum state x e Mn such that p = P(XP% , 
so that 

Tr [«>* (p)] = Tr (p] x Pf) = Tr (p ( x P ( ) = Tr (p) , (51) 

and <&k shows trace preservation on M^. In particular, for k — N, M( = Mn and (51) holds in full generality. 

2) From the same lines of previous point (1) one gets: Tr [<I> (p*)] = Tr (Nfr( ^n/N} = Tr yP^j /N = k/N. 

3) We deduce from Lemma 5 that 

that is positive and Hermitian too. Its normalization, computed in previous point (2), makes this state a sub-normalized 
state. The second statement in (3) is actually the definition of CP unital maps. ■ 

As a particular example, we will consider now the case of fe-extremal CP-TN I maps Ojt associated with rescaled dynamical 
matrices that are product states, cr^ = u <8> for some u e Mn and £ e Vk (&n)- Fot this very special case, the action of 
<£>k on a generic p e Mn reads 

4>(p) = wTr (P ( p) : (53) 

the map sends the entire fc-dimensional subspace spanned by M( into cj, whereas the complementary subspace is 
annihilated. 



7 Concluding Remarks 

In this work we investigated the set of subnormalized quantum states and the set of completely positive, trace non- 
increasing maps. In particular we computed the volume of the set of subnormalized states and provided an algorithm to 
generate such states randomly with respect to the Hilbert-Schmidt (Euclidean) measure. 

We described the structure of the set of CP-TN I maps and computed its volume. It is worth to emphasize that up till 
now no exact result for the Hilbert-Schmidt volume of the set of CP-TP maps is known, (although some estimates can be 
done [29]). On the other hand, in this work we obtained an exact result for the HS volume of the set of TNI maps, which 
includes the set of TP maps. 

Our paper contains several side results worth mentioning. In particular, we found: 

a) the probability distribution of states obtained by partial trace of mixed states of the bi-partite system distributed ac- 
cording to the Hilbert-Schmidt measure; 

b) the volume of the sets of normalized, subnormalized and subtracial states with respect to the measures induced by 
partial trace; 

c) an interpretation of the extremal trace non-increasing maps, which act as trace preserving on certain subspaces. 
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A Proof of Lemma 1 

For any given t e R + , consider the quantity 



r ( \ 2 

G(t,N,K) := © f - X A ' 11 0( A) A '~ N Ft ( Al - A ;) dAl ■ • ■ dAN ' 

J V j=l ) i=\ i<j 



(54) 



where is the Heaviside step function, fulfilling (fx) = (x) for every positive t. By rescaling the variables A, — > tAi 
we get 

r I N \ N 2 
G (t, N, K) = t NK 8 1 - J] Aj Y\ ®W ~" ft _ A J dil ' ' ' diw ' 

^ V ,'=1 J |=1 /</ 

Comparing with (54) gives us the scaling relation G(t,N,K) = t NK G(1,N,K). By taking the derivative at t = 1, we 
obtain 



At 



G (t, N, K) 



NK ■ G(1,N,K) , 



(55) 



and equation (54), together with (9b), yield 



- G (t, N, K) 
at 



N \ N 



[ 6 i ~z a jn® (A '- )A ™n( A '-- A ;) 2 dAi • dA ^ = 7^ 



-N+1,2) 



'=1 " V /=1 / i=l (<; ^-n 

Observe that G ( 1 , N, K) on the r.h.s. of (55) is nothing but the integral in the statement of Lemma 1 , so the result follows. I 



Remark (Volume of M and Al sub with respect to induced measures) 

N N 

The one parameter measure Avk{A\,. . . , A N ) introduced in Lemma 1 is called induced measure [23], 
and appears as natural measure for reduced density matrices of pure states equidistributed on a bipartite 
system of dimension N x K. As it is evident by comparing (7) and (12), the HS-distribution coincides 
with dv N (Ai, . . . , A N ), and with respect to this specific measure we compute the volumes of M N and of 
Al^ ub , given by equation (11) and (15). As a by-product of Lemma 1 , we are able to compute the volumes 
of quantum states and that of sub-normalized states also in the case of induced measure for arbitrary K, 

vowaw^^S! . #( 2,r-»» r<I -" tl,r «-'' (2| - r «- |)r,n 



M 1 MK-N+1,2) 
N 



and 



, x ! Vol[F/fl 

voi V Mf = — yy^ 

V ' N\ NKC ( *~ N+1 ' 2) 



(2n) 



N(N-l)/2 



r(AP) 



T(K-N + l)Y(K-N + 2)---r(K- l)T(K) 



NKT(N 2 ) 
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B Proof of Lemma 2 

Let us define sff := Q>mm = Dim . From (31a) it follows that e (<c) is Hermitian. For any given p e Mn, the hypothesis 



nv mv 



O e CP N m implies p' mm = <D mm p„ v = e { ^ Pnv = TrU^p) < 1 = TrQWp). 

nv x ' 

In other words 

Tr [(l N - e (0) )p] > , M P eM N • (56) 

The matrix (l N - e (4>) ) is Hermitian, thus there exists a unitary U (<s>) such that t/ (0) (Ia, - e (a>) ) t/^ = E ((t) , with H (a>) 
diagonal. Equation (56) must hold for all p e Mn, so that will hold in particular on the sequence of density matrices 

= O '><*• I f/ (0) ' ( 57 > 
Inserting each of the N matrices of (57) into equation (56), and using the cyclic property of the trace, one obtains 

Zf } >Q , V 1 < i < N , (58) 



or equivalently 
Now (32) follows by observing that e (0) = N Tr A o- . 

C Proof of Lemma 3 

Directly from (40-41) we write 



4 0) < 1 , V 1 < i < N . (59) 



f dv K (Ai ,...,A N )= f Y] (1 - A,) ©(A,) Af~ N Y] (A, - A,) 2 dAi . . . dA^ 



By rescaling variables A, — > Ai/N , one obtains the result 

N 



f dv*(Aj , . . . , A N ) = N~ NK f f] (1 - Ad ®{A,) Af~ N \\ (A t - Ajf dA l ... dA N 

J J i=i i<j 

1 r 1 r 1 N ■> 



equivalent to eq. (41) (see [22] for the definition of the Selberg's integral). 



D Proof of Lemma 4 

In Section 4 we decided to express every CP-TNI map O : p i — > p' in terms of a linear superoperator <t> and a related 
dynamical matrix D®, as follows 

P'ik =®ikP j( = D ijPje ■ (60) 

it u 
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According to Lemma 2, the N 2 x N 2 dynamical matrix D® is Hermitian and positive, thus one can benefit of its spec- 
tral decomposition. Defining as |m w ) := Y .. nif. \i; ;') the u th (bi-indexed) eigenvector of Do, corresponding to the 

' I / 'J l J 

eigenvalue m (/i) e R + , the square root decomposition reads 

s s s 

D = 2m w |m w )(m w |=2(V^|m M ))(V^(m w |) = ^|A /i )(A /J | , (61) 

(1=1 £/=l (1=1 

with [aJ. . := Vm^' mff and 5 is given by the number of eigenvalues of different from zero. Equation (61) can be 

L Aij J 

re-expressed in matrix form as 

D 8=ZHM* ■ (62) 

where the symbol x denotes the complex conjugate of x. Now we observe that the each bi-indexed vector [a^,] , of size 
N 2 , can be seen as a square matrix of size N, and reshuffling its indexes one gets [a^J = [Au]^ Equation (60) becomes 

p^= d op^ = ZN p„WL . (63) 



that is the matrix form of (44a). 

The explicit matrix form of the l.h.s. of (44b) reads 



N 

Z K M, - Z Z WLN - Z Z M,M - Z ° r ■ <«> 

//=i /i=i f=i f=i /j=i f=i " 

where we made use of (62), so that (44b) follows from the Proposition of pag. 10. 

To prove the second part of the Lemma, one simply notes that the r.h.s of (61) always defines a positive and Hermitian 
JV 2 x N 2 matrix D, no matter of which S complex N xN matrices have been used. The additional constraint that makes 
D becoming a dynamical matrix, representative of some CP-TNI map O, is given by (34). Imposing (34), in terms of the 
family of matrices A^, is equivalent of imposing (45a), as it can be earned by reversing the chain (64) and by noting that 
& e => o" T e Mff. Finally the claim can be obtained through the same steps (61-63) as before and this ends the 
proof. ■ 
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